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Abstract 

In this paper a technique of compensating for pneumatic 
distortion in pressure sensing devices is developed and veri- 
fied. This compensation allows conventional pressure sens- 
ing technology to obtain improved unsteady pressure mea- 
surements. Pressure distortion caused by frictional attenu- 
ation and pneumatic resonance within the sensing system 
makes obtaining unsteady pressure measurements by con- 
ventional sensors difficult. Most distortion occurs within the 
pneumatic tubing which transmits pressure impulses from 
the aircraft’s surface to the measurement transducer. To 
avoid pneumatic distortion, experiment designers mount the 
pressure sensor at the surface of the aircraft, (called in situ 
mounting). In situ transducers cannot always fit in the avail- 
able space and sometimes pneumatic tubing must be run 
from the aircraft’s surface to the pressure transducer. Re- 
search began at the National Aeronautics and Space Ad- 
ministration Ames Research Center Dryden Flight Research 
Facility, supported by the University of California, Los An- 
geles, to develop a technique to measure unsteady pressure 
data using conventional pressure sensing technology. A 
pneumatic distortion model is reduced to a low-order, state- 
variable model retaining most of the dynamic characteristics 
of the full model. The reduced-order model is coupled with 
results from minimum variance estimation theory to develop 
an algorithm to compensate for the effects of pneumatic dis- 
tortion. Both postflight and real-time algorithms are devel- 
oped and evaluated using simulated and flight data. 

Nomenclature 

Acronyms 

BVP boundary value problem 
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ESP electronically scanned pressure module 

DC direct current 

Symbols 

A c cross-sectional area of pneumatic tubing 

B state equation input geometry matrix 

C observation matrix 

COV covariance 

c sonic velocity 

D pressure tubing cross-sectional diameter 

dx spatial coordinate differential 

G stationary measurement error variance 

G*+ 1 measurement noise covariance matrix 

G{ 1 ) first diagonal component of Gk+ 1 

G( 2 ) second diagonal component of G k + 1 

g k + \ measurement noise function 

I identity matrix 

/ input to deconvolution algorithm 

I deconvolution input estimate 

K *+ 1 Kalman filter gain matrix 

k data frame index 

L length of pneumatic tubing 

N number of data frames available 

O output from pressure sensor 

O deconvolution output estimate 

P a initial local pressure 

i (k prediction error covariance matrix 

P k + x /k+ 1 filter error covariance matrix 

P k + ! f N smoother error covariance matrix 



Pi 

downstream pressure response 

Po 

prescribed input pressure function 

P(x,t) 

pressure function 

Qk* 1 

state noise covariance matrix 

Qc/f, 

prescribed input function error variance 

R 

acoustic resistance of pneumatic 
configuration 

Rk + 1 

residual error vector 

t 

time coordinate 

u k+ 1 

input vector 

V 

enclosed transducer volume 

V' 

effective sensor volume, V + 

X 

spatial coordinate 

X 

state vector 

z 

observation vector 

z 

filter response residual 

Greek 

r fc +i 

prescribed input vector 

Ajt+i 

recursive residual error matrix 


input error vector 

At 

discretization time interval 


dynamic viscosity 

£ 

damping ratio of reduced-order model 

Po 

linearization density 

<J> 

state equation transition matrix 

Wn 

natural radian frequency of reduced-order 
model 


Introduction 

Recent advances in aircraft performance and maneuver 
capability have complicated the problem of flight control 
augmentation. Increasingly, aircraft system designs use 
aerodynamic parameters derived from pneumatic measure- 
ments as control system feedbacks. To be used as control 
feedbacks, pneumatically derived data must be measured 
with accuracy and dynamic fidelity. This has been difficult 
because of pressure distortion caused by frictional attenu- 
ation and pneumatic resonance within the sensing system. 
Typically, most distortion occurs within the pneumatic tub- 
ing used to transmit pressure impulses from the aircraft’s 
surface to the measurement transducer. 

To avoid pneumatic distortion, experiment designers 
mount the pressure sensor at the aircraft’s surface, called 
in situ mounting. In some cases this technique is a viable 


solution. In most cases, as when many pressures must be 
measured in a small surface area or when pressure is sensed 
in an external environment, pressure sensor installations re- 
quire pneumatic tubing from the aircraft’s surface to the 
pressure transducer. 

New non-intrusive pressure sensing technologies, such as 
piezo-film or optical techniques such as laser/Doppler, need 
to be developed. These technologies, however, are still be- 
ing developed and will not be reliably proven for at least 
a decade. Additionally, these technologies are expensive. 
A research effort was initiated at the National Aeronautics 
and Space Administration Ames Research Center Dryden 
Right Research Facility, supported by the University of Cal- 
ifornia, Los Angeles, to develop a technique for measuring 
unsteady pressure data with conventional pressure sensing 
technology. Preliminary results of this research have been 
published in Refs. 1-3. This paper presents further results. 

A pneumatic distortion model, developed in the early 
stages of this research, 2 3 will be reduced to a low-order, 
state-variable model which retains most of the dynamic 
characteristics of the full model. The reduced-order model 
will be coupled with standard results from minimum vari- 
ance estimation theory to develop an algorithm to compen- 
sate for the effects of pneumatic distortion. Both postflight 
and real-time algorithms are developed. Algorithm charac- 
teristics arc evaluated using simulation and flight data. 

Background 

For most flight test applications pressure sensing devices 
are located remotely from the aircraft surface, usually within 
the internal cavity of the aircraft. A connecting length of 
pneumatic tubing transmits pressure changes at the aircraft’s 
surface to the pressure transducer. Pressure variations at the 
surface propagate as waves from the upstream end through 
the pneumatic tubing to the transducer. The wave propoga- 
tion is damped by frictional attenuation along the walls of 
the pneumatic tubing. This damping causes spectral atten- 
uation of the pressure response and produces both a magni- 
tude attenuation and a phase lag. When the wave arrives at 
the downstream (transducer) end of the pneumatic tubing, it 
is reflected back up the tube. Depending on the frequency of 
the incoming waves and the length of the pneumatic tubing, 
the reflected waves either cancel or reinforce incoming pres- 
sure waves. If the waves cancel, further spectral attenuation 
occurs; if the waves reinforce, the power of the incoming 
wave is amplified and acoustic resonance occurs. Spectral 
distortion will not occur in zero volume configurations (with 
in situ mounting) only. This paper will discuss modeling and 
compensating for the effects of internal acoustic dynamics 
of the pressure sensor. External effects, such as orifice dis- 
turbance of the surface flow, will not be modeled. 

A large body of information on the effects of pneumatic 
distortion is available. Early analyses modeled the pneu- 
matic tubing dynamics for specialized inputs 4-9 assuming 
either highly damped or undamped measurement configu- 
rations. While accurately predicting pressure responses for 
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defined types of inputs, these methods could not predict the 
behavior of pneumatic systems subjected to arbitrary pres- 
sure inputs. 

A complex wave model which evaluates the sensor out- 
put for arbitrary inputs was developed in Refs. 2 and 3. The 
model, based on the unsteady Navier-Stokes equations of 
momentum and continuity, was verified by comparisons to 
special case analytical solutions, laboratory data, and flight 
data. The results of this verification indicate that predic- 
tions of the wave model are extremely accurate. Refer- 
ences 1-3 give detailed descriptions of the tests used to ver- 
ify the model. 

Because the wave model is complex, it cannot be easily 
inverted to develop a compensation routine which predicts 
the sensor input (the surface pressure) from the sensor out- 
put (the measured pressure). Thus the wave model must be 
reduced to a simpler, lower order model retaining most of 
the dynamic characteristics of the full wave model. The 
lower order model is more amenable to analytical manipula- 
tion and will be used to derive the compensation algorithms. 

Modeling and Analysis 


problem (BVP), defines the input pressure and solves for 
the downstream response. The basic dynamics equation is 

d 2 P(x,t ) RdP(x.t) 2 a 2 P(i,t) m 

dx 2 

where po is Ihe initial density, P(x, t) is the local pressure 
within the tubing, and c is the local sonic velocity. The pa- 
rameter R is the acoustical resistance which for laminar flow 
is given by 

R = Hfi (2) 

n D 2 

where p is the local dynamic viscosity. The system is as- 
sumed to be at rest initially; thus the defined initial condi- 
tions are 

P(x,0) = P a , and =0 (3) 

ut ji=o 

where P a is an arbitrary local starting pressure and the de- 
fined pressure function at the upstream boundary is 

P«U) = Po(t) (4) 

Satisfying momentum and continuity at the downstream 
boundary 


This section presents the idealized configuration to be an- 
alyzed first. Next, the full wave model of Refs. 2 and 3 will 
be approximated by a simple second-order filter which re- 
tains most of the dynamic characteristics of the complete 
model. The reduced-order model will be written in state- 
variable form and its frequency response will be compared 
with the full wave model’s. 

Idealized Pressure Sensor Geometry and Conditions 

The pressure sensor configuration (Fig. 1) is modeled as 
a straight, cylindrical tube with an axisymmetric volume at- 
tached at the downstream end. The lube is of constant di- 
ameter (D) and length (L). The attached volume (V), rep- 
resenting the internal volume of the pressure transducer, is 
also constant. A longitudinal coordinate (x) starts at the up- 
stream end of the tube, and a time coordinate (t) starts at 
some initial time. 

A small pressure pulse is introduced at the upstream end 
of the pneumatic tubing and propagates down the tube as a 
longitudinal wave. The flow velocity is approximated by the 
radial average, and as a result the pressure within the sensor 
is modeled as a function of x and t only. Since the sensor 
is closed at one end, flow velocities are small and air den- 
sity within the sensor has only small variations. Since the 
flow velocities are small, temperature within the sensor is 
nearly constant. 

Governing Equations 

The governing equations of the full wave model are de- 
rived in Ref. 3. This model, a hyperbolic boundary value 


d 2 P(L,t) . R dP(L,t) t AcC 2 dP{x,t) _ n 

dt 2 + po dt + V dx |,. t 1 

where P(L t t) is the sensor response at the pressure trans- 
ducer. Equations (1) to (5) represent a completely and prop- 
erly posed BVP. The existence and uniqueness of its so- 
lution is demonstrated in Ref. 3. Numerical solutions are 
obtained using implicit differencing in which the pressure 
function is approximated by a series of discrete gridpoints, 
with partial derivative operators approximated by finite dif- 
ferences. The resulting nonsymmetric tridiagonal system is 
solved using traditional elimination techniques. 10 

Reducing the Wave Model to a Second-Order Equation 

Analyses of numerical solutions to the wave model indi- 
cate that for small-amplitude pressure inputs to low-volume 
sensors, the pressure gradient at the downstream boundary 
can be approximated by the spatial average of the gradient 
within the pneumatic tubing without substantially altering 
the dynamic characteristics of the model. If this approxima- 
tion is performed and the BVP is integrated twice with re- 
spect to position, the wave model reduces to a second-order 
equation of the form 3 7 


d 2 Pi(t) 
dt 2 


+ 2 £uj n 


dP L (t) 

dt 


+ UJ 2 P L (t) =^P 0 (O 


( 6 ) 


where oj n is the natural frequency, and £ is the equivalent 
damping ratio. For this model 
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If equation (6) is integrated using a backward zero-order 
hold, 11 the reduced-order model can be written in state- 
variable form as the matrix equation 
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Pl> 

P Lk> I 


0 


1 




1+2£uj,,A t+(u> H At) 2 1 + 2 A f + ( w n A t) 2 

p L>-, \ 

Pik ) 


(u n &t) 2 

l + 2£u>*A t+ (ut»At) z 


p 0k,i 


(7) 


As mentioned earlier, equation (7) retains most of the 
frequency and time response characteristics of the origi- 
nal wave model up through the first wave harmonic. This 
is shown by comparing the analytical frequency responses 
for the full and reduced-order models. Figure 2 shows 
frequency response results for an overdamped sensor with 
pneumatic tubing 8 ft long and 0.06-in. diameter and a negli- 
gible volume transducer. The operating altitude is 40,000 ft. 
The equivalent damping ratio of the reduced-order model is 
1.93, and the natural frequency is 26.75 Hz. Figure 3 shows 
frequency response comparisons for a lightly damped sen- 
sor with the same configuration operating at an altitude of 
20,000 ft. For this configuration, the equivalent damping 
ratio of the reduced-order system is 0.25 and the natural fre- 
quency is 108.93 Hz. 

For both figures the comparisons are good up to the nat- 
ural frequency of the reduced-order model. Beyond this 
both the magnitude and phase angle comparisons begin to 
diverge. The reason for this divergence is clear; the second- 
order approximation does not allow for higher order wave 
harmonics, while the full wave model does. Thus at frequen- 
cies beyond the first harmonic, the reduced-order model may 
not be acceptable and results of the compensation routine in 
the next section will probably not be accurate beyond the 
first wave harmonic. For analysis beyond the first wave 
harmonic a higher order model which allows for additional 
wave harmonics can be used. The model parameters can be 
matched to the full wave model using system identification 
or resonance peak matching techniques. 


Deriving the Compensation Algorithm 


observations of the pressure response at the transducer. In 
general, the technique of infering a system input when given 
knowledge of the system structure and a measurement of the 
system response is called deconvolution. For this applica- 
tion, this is a difficult task. 

If equation (7) is solved for Po i+1 in terms of *v.. 
Pi k , and Pi k ^ y the resulting equation is numerically ill- 
conditioned; small changes in P Lk produce large changes 
in Po^ ■ As a result, noise in the pressure measurements at 
the transducer end of the sensor is over amplified by the nu- 
merical instability of the computation and overwhelms any 
estimate of the input pressure. The measurement noise must 
be identified and controlled as a part of an estimation loop. 
The algorithm that performs the error control function, pre- 
sented in this section, will be used to infer the external pres- 
sure input to the system based on observations of the pres- 
sure response at the transducer. 

First, a general deconvolution algorithm as developed in 
Ref. 3 will be presented. The technique, developed using the 
techniques of minimum variance estimation theory, is simi- 
lar to results presented in Refs. 12-14. Both postflight and 
real-time versions of the algorithm will be presented. After 
the general algorithms are presented, they will be applied to 
the pressure compensation problem. 

Minimum Variance Deconvolution Algorithm 

The deconvolution algorithm defines an input to the state- 
variable model and identifies the residual between the pre- 
scribed input and the actual input using minimum variance 
estimation techniques. All notational conventions in this pa- 
per follow the conventions of Ref. 15. The assumed form of 
the model is 


x k +i = d >x* + B[r* +1 + fi k+1 ] (8) 

= CX jt +1 + g k+ 1 (9) 

where, 

Uk+ i = r* +1 + <5fc+i (10) 

In equation (10) [/*+ j is the actual input to the system, r*+i 
is the defined (nonrandom) input to the system model, and 
& k + 1 is a random process which models the residual between 
the actual and prescribed inputs. The covariance of 8 k+ \ is 
defined as 

COV[5 fc+l ] = Q* +1 

and is assumed to be known. 

Measurements of the state vector ( Z k+ 1 ) are observed in 
the presence of noise. The measurement noise covariance is 
given as 

COV[g k+l ] = O k+l 
and is assumed to be known. 


The compensation algorithm will be developed in this 
section. The objective is to develop an algorithm to infer the 
external pressure input to the measurement system based on 


For notational convenience, the state transition matrix 
(d>) and the input and measurement geometry matrixes ( B 
and C) are written as constants. 
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Two complete estimation algorithms will be presented. 
The first, intended for postil ighl data reconstruction, is im- 
plemented as a two-pass forward filler, backward smoother. 
As developed in Ref. 3, the postllight smoothing algorithm 
relies on the use of all available data measurements, both 
past and future. The algorithm, a hybrid of the Kalman filter 
and the Rauch smoother, 15 uses information in the Kalman 
filter innovations to identify differences between the defined 
input and the true input during the reverse pass of the algo- 
rithm. The second algorithm, intended for real-time opera- 
tion, will be derived from the smoothing algorithm and im- 
plemented as a time-recursive filter. This algorithm docs 
not use information from future observations of the sys- 
tem response. 

Postflight Smoothing Algorithm. Assuming that ob- 
servations of the system response arc available for data 
frames k = 0 , 1 , . . . , N, the following sequence of com- 
putations are required for the smoothing algorithm: 

Kalman filter loop, performed for A: = 0 , 1 , . . . , N 

X k+ 1/* = QXit/ic + BFVi 

P k +i Ik = ^Pk/k® 1 + BQk+] B l 

Z k +\ /k = Z k + 1 — CX k +\j k 
K k . , = Pt + i/tC T [CP*.i/*C T + G fc+1 r l 
= X k+ \/k + Kk+\Zk+\jk 

P fc+ 1 /fc+ 1 = f I ~ Kk+ 1 C ] Pk+\/k 

where P k +]/k is predicted filter error covariance, 
^fe+i/A+i is the adjusted filter error covariance based on the 
current observation, K k +\ is the Kalman filter gain matrix, 
and Z k+ \/ k is the filter innovation. The matrixes BQ k +\B l 
and 1 arc the assumed covariances for the input and mea- 
surement noise processes. The vectors X k + \ j k and X k + \ j k + i 
are the predicted and corrected estimates of the filter state 
vector. Resulting estimates for the stale vector, the error co- 
variances, and the innovations process must be stored for 
each recursion, k = 0 , 1 , . . . , N. 

At this point the input residual estimation loop is per- 
formed backwards for A: = Af, M — I , . .. ,0, using 

stored data 

Pk/N = A k Rk+\/N + O l \CP k j k -\G x + G k 1 x Z k f jt_i 
&k/N ~ QkB T R k fN 

where R k /u is the deconvolution residual vector, defined by 
Pk/k - 1 Rk/N ~ X k /tf - X k / k _\ ( 1 1 ) 

and 

A* 1 ' = ( / - K k C)' r <t> r = I <t>U - K k C) ] T 
Finally, the input estimate is computed as 
UkjN ~ PkjN + &k/N 

The final value for the residual vector, R n +\ini is given by 
equation (11) with k replaced by N 

Rn+[/N = I Pn+\/n] 1 l/N ~~ X-M+\ftf] - 0 


Real-Time Filtering Algorithm. The real-time algo- 
rithm relies only on past and present data to perform the 
estimation. Assuming that observations of the system re- 
sponse are available for data frames k - 0 , 1 , . . . , fc + 1 , the 
residual estimation is performed with k + 1 replacing N 

h+\/k+\ = Qk+\B T R k+ \/ k +\ 

i’U+l/A+l = — X k +]/ k ] 

therefore, 

^fc+l/Jfc+l = Qk+\ B^ [ P k + 1 j k ] f^fc+l/fc+1 1 /fcl 

and substituting in the Kalman filter 

i /k + 1 = Q/c+i B T [ P k +\/ k ] 1 K k +xZ*\,k 
Applying the definition of the Kalman gain matrix 

S t+ . /*. 1 = Q 1 B T C T [ C P k+ 1 C T ] - 1 z k + , /* 

The input may be estimated as a part of the Kalman filtering 
loop according to the sequence of computations: 

Kalman filter step performed for fc = 0 , 1 N 

-X’t+l/i = 'bXk/k + STi.i 
Pk,Uk = ^Pk/k^ T + SQit+i S T 

Z k +\j k = Z\ k+i — CX k +\j k 
K k , , = P k + ]/k C J {CP k+ i fk C r + G* +1 ]-’ 

X k + 1 /fc+ 1 = X* +1/ * + K k+ \Z k +\/ k 
P k + 1/Jk+l “ [ f “ Kk+\ C]P k +\/k 

Residual estimation step 

= Q k +\B T C T [CP k+ iC T ]- l Z ki . ]/k 

and 

U k +\/ k +\ = Pk+\/k+\ + <5*+i/fc+l 

Applying Minimum Variance Deconvolution 
Algorithms to the Pressure Compensation Problem 

The above algorithms will now be applied to the pneu- 
matic distortion problem. Equation (7) may be written as 
the first-order matrix difference equation 

X k + 1 =<t>X k + BU k + 1 (12) 

If the downstream pressure is observed in the presence of 
noise, then 

Z k + 1 = X k +\ + g k +\ ( 13) 

If the substitution 

Ui t+i = r k +\ + Sjt+i 
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is performed, equations (12) and (13) are the same form as 
the assumed deconvolution model, equations (8) and (9), 
and the deconvolution algorithms may be directly applied. 

The input function (T k+X ) is determined by solving equa- 
tion (7) for the input pressure in terms of the downstream 
pressures, and evaluating the result using state estimates 
from previous time recursions of the filter loop. Amplified 
noise in T * + 1 is estimated and controlled by the deconvolu- 
tion algorithm through judicious selection of the state and 
measurement noise covariance matrixes. 


The state noise covariance matrix is derived from equa- 
tion (7) as 


BQ k+ ,B r 


Q rr 

WUk*i l + 2£w*At+(u/,,At) 2 


(14) 


where Qy k , x is the variance in the residual <5^+1 . Similarly, 
assuming that the measurement noise is white, the measure- 
ment noise covariance matrix is 


G *♦! 


G( l) fc+ i 0 
0 G(2) fc+1 


(15) 


where G{ l)*+i is the error variance of the measured down- 
stream pressure at data frame k and G( 2)* + i is the error 
variance in the measured downstream pressure at data frame 
it + 1 . For stationary measurement processes 


G( l)* + i = G( 2)* +1 = G 


Results of Deconvolution Analysis 

Selected results of the minimum variance deconvolution 
analysis will now be presented. Filtering tuning will be 
demonstrated using simulated data. Next the deconvolution 
algorithms will be applied to flight data. Both the postflight 
smoothing and real-time filtering algorithms will be ana- 
lyzed and compared with reference pressure data obtained 
from flight experiments. 

Filter Tuning 

Selecting proper assumed values of the error covariance 
matrixes is critical. This practice, referred to as “filter tun- 
ing,” is typical of minimum variance estimators. The appro- 
priate values for the covariance matrixes are usually unique 
to each configuration. Proper tuning of the deconvolution 
algorithms will be shown using simulated data generated 
from numerical solutions of the full wave model. Simulated 
data is used because the actual measurement and input er- 
ror covariances can be predefined. The simulated data were 
generated using the full wave model assuming the configu- 
ration associated with the overdamped case of Fig. 2; a sen- 
sor configuration with pneumatic tubing 8 ft long, and of 
0.06-in. diameter, a transducer with negligible volume, and 


operating at an altitude of 40,000 ft. For this configuration 
and altitude, the equivalent damping ratio of the reduced- 
order system is 1.93; the natural frequency is 26.75 Hz. Sim- 
ulated zero-mean white noise with a variance of 0.1 ( ^r) 2 
was superimposed on the output time history; simulated 
zero-mean white noise with a variance of 1 .0 ( Tpr) 2 was 
superimposed on the input time history. The input and out- 
put noise were uncorrelated. 

The minimum variance deconvolution analysis was sim- 
ulated using the postflight smoothing algorithm. Figure 4 
shows the simulation process. A defined input is superim- 
posed with additive noise and the result (I) is input to the 
wave model. The output from the wave model is superim- 
posed with additive noise to give a pseudomeasurement (O) 
time history. Using values of the natural frequency and the 
equivalent damping ratio provided by the wave model and 
prescribed values for the input and observation covariances, 
the pseudomeasurement is used to perform the deconvolu- 
tion. The resulting deconvolution estimates of the model 
input and response are indicated by the symbols I and O. 

Figures 5 and 6 show the effects of selected values of G 
on the estimates of U. In these figures, the value of Q y was 
fixed at the appropriate values of 1.0, while the value for G 
was varied. Recall that the actual measurement error vari- 
ance is 0 .1 ( J£r) 2 . Figure 5 shows the deconvolution re- 
sults for G = 10 , a value which is much too large. Three 
time history graphs are shown. Time history overplots of 
the measured downstream pressure response (O) and the 
smoother estimate of the input pressure (/) are shown in the 
bottom graph. Time history overplots of O and the estimated 
response (d) are shown in the middle graph, and time his- 
tory overplots of the actual input (I) and I are shown in the 
top graph. Comparisons are poor. This set of weights does 
not allow the estimation algorithm to extract enough infor- 
mation from the measured pressure data; resulting estimates 
are overly smoothed. 

The above situation cannot be improved by arbitrarily 
lowering the value for G . Figure 6 shows similar time his- 
tory comparisons for a value of G which is much too small, 
0.01. The end result is a virtually perfect match between 
the measured and predicted response. As discussed earlier, 
this results in over-amplification of the measurement noise 
and the estimate is overwhelmed. The selected value for G 
did not allow the deconvolution algorithm to identify and 
control noise in the measured response. The filter is essen- 
tially trying to perform an open-loop deconvolution of the 
state equation. 

Figure 7 shows the effect of selected values of Qy on the 
estimation of U . The value of G was fixed at the appropriate 
value of 0.10, while the value for Qy was set at 10. The ac- 
tual value for Qy is 1.0. The resulting estimate contains ex- 
traneous harmonics that are not present in the original input 
signal caused by overamplification of mid-range frequen- 
cies. The dynamic constraints of the reduced-order model 
have not been enforced strongly enough. 
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The appropriate level of smoothing and compensation is 
shown in Fig. 8. Here the proper values for the measure- 
ment and input error covariances arc used. In this ease the 
estimated input matches the actual input well, and much of 
the superimposed noise has been removed. 

Application to Flight Data 

The deconvolution algorithms will now be applied to 
flight data from the experiments described in Refs. 1-3. The 
flight test configuration mounted on the wing of a high- 
performance aircraft (Fig. 9) consisted of a set of static lest 
orficies connected by pneumatic tubing to a single electron- 
ically scanned pressure module (ESP), and a single in situ 
mounted pressure sensor. The test matrix had ports of 0.02, 
0.04, and 0.06 in. in diameter and allowed pneumatic tub- 
ing sections to be interchangeable. Reference pressure data 
were provided by a very low volume piezioelectric wafer 
pressure transducer glued to the aircraft skin adjacent to the 
test matrix. Since the in situ transducer was mounted on the 
aircraft skin, it was exposed to temperatures considerably 
lower than the recommended normal operating range. This 
caused the direct current (DC) output to drift off the scale. 
To eliminate the drift, the reference sensor output was elec- 
tronically coupled with a third-order band-pass filter (1.16 
to 101.7 Hz) to remove the DC component. The DC com- 
ponent was restored to the reference data using the output 
from a pressure port connected to the ESP module by 6 in. 
of pneumatic tubing. The output from this port was low- 
pass filtered using a second-order digital filter with a rolloff 
frequency at 1.0 Hz. The two filtered signals were merged 
using a complementary filter. 

Pressure data obtained from the test matrix measurements 
were compensated using both the postflight smoothing and 
real-time filtering algorithms. The results will be compared 
to the reference data reconstructed using the complemen- 
tary filter. The flight data deconvolution process (Fig. 10) 
is similar to the simulation process described earlier. Two 
flight data cases will be presented; one for an overdamped 
configuration, and one for an underdamped configuration. 

Deconvolution results for the overdamped flight maneu- 
ver obtained at 42,000 ft altitude with a sensor configuration 
with 8 ft of 0.06-in. -diameter pneumatic tubing and negligi- 
ble transducer volume are shown in Fig. 1 1 for the smooth- 
ing algorithm and Fig. 12 for the filtering algorithm. For this 
configuration and altitude, the equivalent damping ratio of 
the reduced-order system is 2.125; the natural frequency is 
26.75 Hz. Again, three time history plots arc shown. The 
lower graph compares 7 with O. The middle graph com- 
pares O and 6. The upper graph compares I to /. The 
comparisons arc excellent for both the postflight smoothing 
and real-time filtering algorithms. The estimated input does 
not contain the high-frequency noise present in the refer- 
ence data. 

Figure 13 shows deconvolution results for an under- 
damped flight maneuver at 24,000 ft altitude with a sen- 


sor configuration having 2 ft of 0.06-in.-diameter pneumatic 
tubing and negligible transducer volume for the smooth- 
ing algorithm. Figure 14 shows results for the filtering al- 
gorithm for the same configuration and altitude. For this 
configuration and altitude, the equivalent damping ratio of 
the reduced-order system is 0.279; the natural frequency is 
107.18 Hz. Both the smoothing and filtering algorithms 
work well. The smoothing algorithm removes nearly all 
of the resonant harmonic, approximately 100 Hz, from the 
measured pressure response. The filtering algorithm re- 
moves a sizable portion of the resonance. Both algorithms 
account for phase lag well. 

These results are typical of most flight maneuvers ana- 
lyzed. Both the postflight smoothing and real-time filtering 
algorithms performed well. All tuning of the filter covari- 
ances was perform off-line using preflight information and 
arbitrary selection. No on-line adaptation of the covariances 
was performed. 

Summary and Concluding Remarks 

In this paper, a general numerical technique for obtaining 
unsteady pressure measurements using conventional pres- 
sure sensing technology was developed. A mathematical 
model, based on the Navicr-Stokes equations of momen- 
tum and continuity, was reduced to a simple state-variable 
model. Frequency responses of the full and reduced-order 
models were shown. 

Using a collection of standard results from minimum vari- 
ance estimation theory, the reduced-order model was incor- 
porated into a deconvolution algorithm which estimates the 
system input based on the system output. Both postflight 
smoothing and real-time filtering algorithms were devel- 
oped. Overamplification of measurement error is controlled 
by careful selection of the elements of the algorithm state 
and measurement error covariance matrixes. 

Performance characteristics were discussed. Filter tun- 
ing characterises were demonstrated using simulated data. 
Applications of the techniques to flight data were presented. 
Both the postflight smoothing and real-time filtering algo- 
rithms performed well for overdamped and underdamped 
sensor configurations. 

The results of this research are promising for the model- 
ing and compensation techniques developed. These tech- 
niques provide a means of predicting and accounting for 
the effects of pneumatic distortion in pressure sensing de- 
vices. They also offer a reliable and cost effective means of 
measuring unsteady pressure data. These techniques will be 
useful until new pressure sensing technologies can be devel- 
oped and proven. 
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Fig. 4 Deconvolution simulation scheme. 
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Fig. 5 Effect of improperly selected measurement error covariance, Qy = 1 .0 , G = 10 .0 . 
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Fig. 7 Effect of improperly selected input error covariance, Qu = 10 .0 , G = 0 .1 . 
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Fig. 8 Effect of properly selected input and measurement error covariance, Qu = 1.0, 
£7 = 0.1. 
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Fig. 10 Flight data deconvolution algorithm. 
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Fig. 1 1 Flight maneuver time history comparison for the postllight deconvolution algorithm 
output for an overdamped system. 
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Fig. 12 Flight maneuver time history comparison for the real-time deconvolution algorithm 
output for an overdamped system. 
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Fig. 13 Flight maneuver time history comparison for the postflight deconvolution algorithm 
for an underdamped system. 


17 





Time, sec 


Fig. 14 Flight maneuver time history comparison for the real-time deconvolution algorithm 
for an underdamped system. 
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